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Abstract. A massive black hole at the center of a dense stellar system, such as a globular cluster 
or a galactic nucleus, is subject to a random walk due gravitational encounters with nearby stars. 
It behaves as a Brownian particle, since it is much more massive than the surrounding stars and 
moves much more slowly than they do. If the distribution function for the stellar velocities is 
Maxwellian, there is a exact equipartition of kinetic energy between the black hole and the stars 
in the stationary state. However, if the distribution function deviates from a Maxwellian form, 

■ the strict equipartition cannot be achieved. The deviation from equipartition is quantified in 
this work by applying the Tsallis q-distribution for the stellar velocities in a q-isothermal stellar 

, system and in a generalized King model. 
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(T) , 1. Introduction 

A massive black hole as a Brownian particle conducts to a equipartition of kinetic en- 
ergy at the equilibrium, when the distribution function of stellar velocities is Maxwellian. 
When it is non-Maxwellian is not surprising to find a deviation from equipartition. This 
deviation is defined as r\ = (M (V 2 )) / (m(v 2 )) , where M and V represent the mass and ve- 
locity of the black hole, and m and v masses and velocities of stars. Following Chatterjee 
et al. (2002) we have r\ = [3 J™ /(r, v)vdv J* °° f(r, v)v 2 dv}/[f{r, 0) / °° / (r, v)v A dv], with 
f(r, v) representing the distribution function for stars field. Evidently, for Maxwellian 
distribution rj = 1. However, astrophysical systems are non extensive since they are 
subject to long range interactions. Boltzmann-Gibbs (BG) thermostatistics is able to 
predictions in extensive systems, in the sence that microscopies interactions are short 
or ignored, temporal or spatial memory effects are short range or does not exist and is 
valid ergodicity in the phase space. BG domain is enlarged by non extensive statistical 
mechanics, based on Tsallis entropy (Tsallis 1988) and in the g-distribution function: 
p(x) oc [1 — (1 — q)/3x] 1 /( 1 ~ q \ Boltzmann distribution is recovered in the q — > 1 limit, and 
all the usual statistical mechanics as well. A summary about mathematical properties of 
these functions can be found in Umarov et aZ.(2006). 
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2. Extended stellar models, deviation from equipartition and results 

The isothermal sphere is the simplest model for spherical systems. It was generalized by 

Lima & Souza (2005) based on the distribution function f q (v) oc [l — (1 — q)v 2 / (2a 2 )] 1-9 
and the deviation results in r/(q) = (7— 5g)/[2(2 — q)]. Most commonly used are King mod- 
els, based on truncated isothermal spheres. Extended King models were presented in Fa & 
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Figure 1. At the left, deviation for q-King model with q = 0.94, and q = 0.8. The limit value 
is rj = (7 — 5g)/2(2 — q). At the right, r\ for q = 1.05. King means 5=1, the usual King model. 
Strict equipartition implies rj = 1. 

Pedron (2001) and applied to fit surface brightness of the NGC3379 and 47TUC with ex- 
cellent results. The distribution function in g-King models is fk q (e) = Pi(27rer 2 )~ 3 / 2 {[1 + 
(1 — q)e/t7 2 ] 1 /( 1 ~ q ) — 1} for e > and fk (e) = for e ^ 0. Here e is the relative energy 
e = ip — v 2 /2 > and tp = —(f>(r) + 4>o the relative potential. The central potential is 
W = ip(0)/a 2 and in the 4>(0)/a 2 — > oo limit the isothermal sphere is recovered. For such 
model the r\ value for q < 1 and w — —<j>(r) is 



where (3 is the incomplete Beta function, a 
q > 1 wc obtain 
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and A = 1 



(2.1) 
aw. For 

(2.2) 



where b = q — 1, A = 1 — bw, and a;m = x ^_ w ■ Furthermore, ii = f Q m (l + x) 1 / b dx, 

h = !a m {l + x)- l ' b x l ' 2 dx : and J 3 = /""(l + x)-V b a?/ 2 dx with the conditions A = 
l-bw^{) and o < 1 + - . 

In the Figure [1] results indicate that there is a deviation from equipartition a priori, 
even at very long time-scales. The equipartition is never achieved in both cases ( q < 1 
and q > 1). The g parameter will be, in some way, dictated by the system. 
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